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We introduce a new method to detect the absolute neutrino mass scale. It uses a macroscopic 
mass of tritium source. We explain that the neutrino mass can be measured by scaling the mass 
difference of the source between initial and final state, and its heat value. This method is free from 
the electron energy resolution limit and the statistical error. We estimate the required accuracy to 
measure the neutrino mass. We also report that the {u, d} component of the CKM matrix, \Vud\ 
may be determined in 10 -6 accuracy as an application of this work. 
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I. INTRODUCTION 



The neutrino was introduced by E. Fermi in 1935 [lj. Since then, the neutrino is treated as the massless fermion. 
However, the neutrino oscillation demands at least two of three neutrinos to have non-zero masses, and the squared 
mass splittings are 0] 

' Am 2 = (6 - 9) x 1(T 5 [eV 2 ]. 

- , (1) 

I ; Am 2 3 = (1 - 3) x 1(T 3 [eV 2 ]. 

Ph| 

Meanwhile, the neutrino mass has been searched by many experimentists. Today, one of the most established 
i ' method is as follows: Prepare the tritium (T) sample. Tritium decays into helium-3 ion, electron, and anti-electron- 
neutrino. This process is written as 

m ■ 

>' T -> 3 He+ + e~ + D e . (2) 

m ; 

Then detect the electron energy distribution. The endpoint spectrum of this distribution depends on the neutrino 
mass. So, you can determine the neutrino mass by detecting the endpoint, accurately. The latest upper bound 
is m v < 2 eV [2,]. Some other methods are performed, for example, the neutrinoless double beta decay and the 
cosmological observation, and they give the upper bounds {ra v ) < 0.7 — 2.8 eV Q and Y,m u < 2.0 eV Q, respectively. 

If the mass hierarchy is normal, the lower bound of second heavy neutrino mass is > 0.008 eV as the lightest 
neutrino is massless. If we have this sensitivity, we surely determine the absolute neutrino mass. 

However, the KATRIN experiment, which will start in 2012 Q, is designed with sensitivity to measure the effective 
neutrino mass m v > 0.2 eV. 

k> , The absolute neutrino mass scale is one of the undetermined parameter of the Standard Model (SM). Some models 
beyond the SM, for example, GUTs [(| 0] predict the absolute neutrino mass scale. Moreover, for reliability, the 
alternative method to detect the neutrino mass is important. 

Our new method explained in this paper is one of the tritium beta decay experiments. However, today's tritium 
beta decay experiments are suffered from some difficulties as follows: First, the required electron energy resolution is 
more and more severe. Next, the produced electron energy spectrum is distorted by random multiple scattering in 
the source. These difficulties are caused by detecting the electron kinetic energy. 

The KATRIN experiment [5J uses the huge detector which spectrometer is 23 m long and the gaseous tritium source 
will consist of a 10 m long. Our method needs less space since we measure the mass difference and heat value of 
macroscopic sample as explained later. 

This paper is organized as follows: In Section UH we explain the new method. In Section [ill! we derive the required 
accuracy to measure the neutrino mass. In Section [V] we discuss the results and summarize this work. 
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FIG. 1: The basic concept of the new method. We detect the mass deference and the heat value. 



II. THE NEW METHOD 



The concept of the new method we introduce is as follows: Elementally, the beta decay process is written as 
?i° — > p + + e~ + v e . In this process, as explained later, the v e mean energy (E v ) depends on the effective D e mass m„, 
which is explained in Appendix |A"1 Then, we can determine ra v by detecting (E v ). 

This method has some advantages compared to the existing ones. In final state of beta decay, v e has no interaction 
with the sample. Since {E v ) is the mean value, we don't have to measure the beta decays event-by-event. Then, we 
can enlarge the mass of the sample to the macroscopic scale. 

We employ the tritium source because the half-life of tritium is approximately 12.33 years; energy difference between 
the initial and final state particles is small; tritium is easy to obtain; tritium has a small nucleon number. The long 
lifetime enables us to detect the mass difference and the heat value precisely. The small energy difference suppresses 
the increasing temperature (as explained in Appendix|B]), and all of the produced electrons are captured in the sample. 
Small nucleon number corresponds to the large events par unit mass. 

To determine (E„) experimentally, we introduce YiE u as the total energy of produced neutrinos during the exper- 
iment, and N as the number of beta decay events during it. Using them, we have {E v ) = TiE u /N . Since we cannot 
detect YiE v directly, it must be treated as a missing energy. Then, we consider the energy conservation between the 
initial and final state sources. In the initial state, we set the condensed tritium which have the mass Mj in a container. 
It includes some impurities. The mass of the container and the impurities is represented as %. In the final state, 
some of tritium decay into 3 He + + e~ + v e . v e have no interaction and go away. 3 He + and e~ have too small energy 
to penetrate the container (typically they have some keV momenta), then they stay in the container and construct 
the 3 He atoms. The sum of 3 He produced in the experiment and the remained tritium have mass Mf. The sample 
releases thermal energy Mq since the beta decay is an exothermal reaction. Assuming that the sample is isothermal 
during the experiment, Mq is dissipated away from the sample by the heat conduction and radiation. Therefore, the 
total mass in the final state is Mf + x- Then, the energy conservation between the initial and final states gives 

Mi + x = Mf+x + M Q +T,E l/ . (3) 

On the other hand, the number of decay event N is given by 

AM 

N= , (4) 

tot - torc 

where AM = Mj — Mf, and tot and TOn e are the masses of tritium and helium-3 atoms, respectively. Using these 
equations, the mean neutrino energy is given by 

VE V , , f, Mq 



(E v ) = ^ = (mr- m Hc ) ( 1 " ±^ ) ■ (5) 

Therefore, we have to detect tot — "t-Hc, Mq, and AM to determine (E v ). First, tot — tojjo should be given by 
other experiments as explained in Appendix |D] Second, to detect Mq, we introduce the technology of differential 
scanning calorimeter ||. Fig. [T]is a schematic drawing of the new method. S is a test sample and R is a reference 
one. S and R are located symmetrically. R releases the thermal energy controlled by the electric heater. TM is the 
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thermo-module and TS is the thermo sensor which measures the difference of released thermal energy between S and 
R. Measuring the electric energy which is injected to cancel the temperature difference between S and R, we can 
determine the Mq as released thermal energy from S. Last, AM is detected by the weighting machine. 

Now we study the neutrino mass dependence of (E u ). The leading order differential decay width which contains 
exact m v effect takes the form 



dT _ G 2 F \V u d\ 2 y(g ~ ml){(M 3 - 2m T £„) 2 - 4mg(M 4 - 2m T E u )} 
'd~E u ~ (2tt) 3 M a - 2m T E v 

x ((1 + C A fE v {M x - 2m T E v ) + (1 - Ci)m He+ (M 3 - 2m T E u )\ 

L l J (6) 

+ 5 {(1-CU) M 2 - 2(1 -C A )m He+ m T j ^ _ 2mT ^ 

_ n _ r x 2 f 3(A/3-2m T ^) 2 (m T ^g,) 2 (E 2 : - ml){(M 3 - 2m T E„) 2 - 4m 2 e (M 4 - 2m T E„)} } 
( C A J m T | 6(M 4 -2m T ^) 2 6(M 4 - 2m T E u ) 2 J 

where Gf, V^d, Cyi, m e , and m He + are the Fermi constant, the {u,d} component of the CKM-Matrix, the axial 
current coupling constant, electron mass, and helium-3 ion mass, respectively; and 



The neutrino mean energy is given by 
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(7) 



dE, 



d£„ = / dE v , (9) 



E Max = Mi - 2m e w He + , 1Q s 
2tot 



III. THE REQUIRED ACCURACY 

is affected not only by m v but also by Ca, m Vl m Ho +, m e , and tot — to Hc + — TO e- How accurate must we 
measure these quantities to detect the neutrino mass? To estimate this, we introduce 

jj /e , / xr y x x sv (E„(m v ;C A + SC A , to Ho + +<5m Ho +,m e + <5m e )} - (-E„(0; Ca, m Ho+ ,rn e )} 

« ro„;iC 4) imHe+,dro e = — — ; - — - TF< — = — — rr. (11) 

(fi„(m„;GA + d6A,m He + +om Ho +,m e + om e )} + {E v {0; Ca, m He +, m e )} 

Here we use the quantities, Ca — 12.6, to Hc + = 2800 MeV, m e = 0.511 eV, tot — m Hc+ — m e = 0.2 MeV, and 
5Ca, Smn c + , Sm e mean the errors of Ca, toh + , m e, respectively. 

Actually, Ca, m^ c +, m e , and tot — to Hc + — rn e given above are not the true values, and also the Eq. does not 
give the accurate decay width since it contains no corrections. However, it does not matter for our purpose. We just 
want to estimate the order of m v effect and errors caused by SCa, <5to-Hc+! and so on. 
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We first consider {E v ) differences caused by m u as 

5(E U (0.1 eV; 0,0,0)) ~ 1.12 x 10 -11 , 
5{E V (1 eV; 0, 0, 0)) ~ 1.12 x 10~ 9 . ^ 

The required accuracy on (E v ) increase by two digits as the neutrino mass which we want to measure becomes 
one-tenth. 

On the other hand, the normalized error caused by Ca is estimated as 

5{E V (Q; 10~ 4 , 0, 0)} ~ 2.44 x 10" 11 , 

S(E v (0; 10~ 2 , 0, 0)} ~ 2.43 x 10~ 9 . ^ 

The required accuracy on Ca increases by one digit as the required accuracy on (E v ) increases by one digit. The 
reason why the required accuracy on Ca is milder than that on (E v ) is explained in Appendix [Cl 

After all, we have to measure Ca in 10~ 4 (10~ 2 ) order to detect 0.1 eV (1 eV) neutrino mass because to detect the 
neutrino mass excess, the Ca originated error in (E v ) has to be less than the (E v ) variance caused by m v . 

A. Other Uncertainties 

To measure m„ accurately, we must overcome some other uncertainties. They are namely mu e +, m e , tot — TOjj e + ~ 
m e , Mi - (Mf +Mq)- E£ 7 , and M Q + HE 1 . 

First, we consider 6(E V (0; 0, 5m Ho +, 0)). Fixing other variables, it is given by 

5{E U (0; 0,2800 x 10" 8 MeV,0)} ~ -1.54 x 10~ 14 , 

S(E v (0; 0, 2800 x 10~ 10 MeV, 0)} ~ -1.54 x 10~ 16 . ^ 
Next, we consider 5{E U (0; 0, 0, 6m e )). Fixing other variables, it is given by 

5(E U (0; 0,0, 0.511 x 10~ 8 MeV)) ~ -1.45 x 10" u , 

6{E U (0; 0, 0, 0.511 x 10" 10 MeV)) ~ -1.45 x 10~ 13 . ^ 5 " > 

Third, we consider (E v )' = {E v {rriT — to Hc + — to 6 ))/(tot — to Ho + —m e ) error caused by Tot— to Hc + —m e uncertainty. 
Fixing other variables as before, it is given by 

(£„(0.02 x (1 + 10- 8 ) MeV))' - (£„(0.02 MeV))' u 

(^(0.02 x(l + 10- 8 ) MeV))' + (^(0.02 MeV))' ' ' 

(E v (0.02 x (1 + IP' 10 ) MeV))' - (^(0.02 MeV))' 13 

(^(0.02 x (1 + 10- 10 ) MeV))' + (^(0.02 MeV))' 

It is reasonable to consider (E u )' because in Eq. ([S]), tot — toh varies along with (E v {m^ — to Hc + — Tn e )), and 
tot — TOHe is easy to calculate from tot — ""Ihc+ — m e- 

These errors are proportional to (E v ) error. If we want to measure m v to 0.2 eV as the KATRIN goal, (E v ) has to 
be determined with 5{E V ) < 10~ n , and then <5to Hc + /to Ho + < 10 -5 , Sm e /m e < 10~ 8 , and 5(mr — to Ho + — m e )/(mT — 
to Hs + — TO e ) < 10 -8 , respectively. 

Last, AM and Mq in Eq. (|S|) require the same accuracy as (E v ). 

As a result, the relations between m v which we want to detect and the required accuracies of related variables are 
shown in Table fl] 

The required accuracy in Ca is comparatively low. This is because Ca in Eq. © does not appear in the leading 
order of E~~ ax . This situation is true for other corrections, for example, radiative correction and nucleus-dependent 
correction. The errors in the Fermi constant and CKM- matrix do not affect (E v ) at all. The Fermi function for the 
Coulomb correction is calculable. 



IV. \V ud \ MEASURMENTS 



As an application of this work, we can determine T (lifetime), and Ca accurately up to theoretical uncertainty. 
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TABLE I: The relation between m v which we want to detect and the required accuracies of related variables. For example, if 
m„ = 0.2 eV, then we have to determine Ca as 5Ca/Ca < 10~ 4 , also m He + as 5m He +/m Ho+ < 10~ 5 , etc. 



m v we want to detect 


20 eV 


2 eV 


0.2 eV 


0.02 eV 


0.002 eV 


{E v ) 


io- 7 


1Q -9 


io- 11 


ht 13 


nr ls 


C A 


1 


10~ 2 


10~ 4 


1Q -6 


io- 8 


m Hc+ 


io- 1 


IO -3 


10" 5 


IO" 7 


10" 9 


m e , m T - m He+ - m e 


io~ 4 


10~ 6 


io- 8 


io- 10 


1Q -12 


AM, M Q 


IO" 7 


IO" 9 


1Q -ii 


1Q -13 


10" 15 



(17) 



The number of decayed tritium N\ and at time t\ and t 2 , respectively, are given by 

m = iv (i- e - rti ) 

N 2 =N (l-e- rt *), 

where iVo is the number of tritium in the initial state. Then, defining the ratio N2/N1 = A, T and A are related as 

dA ( rt ie - rt i _ Tt 2 er Tt i \ dF 

'- rtl 1 -^ iJr (18) 



.4 



1 



1 



r 



- ^(h-h)—, 

where we use the approximations Tt\ -C f and Tt2 <C f in the second line. Therefore, if we measure dA/A = 0(10 ) 
for a two-year experiment, the T ambiguity is dT/T — ©(lO^ 6 ) since T(t2 — ti)/2 ~ 10 -1 . We note here that A is 
actually determined without mx — tihc value since A can be written as 



A = ^ 



Mi-M f2 

771 T — 771 H c 



AM 2 
AMi' 



(19) 



where the indices 1 and 2 represent the values at time t\ and t%, respectively. 

Ts contain \V u d\ and Ca as parameters. If we can determine C'a with 10 -6 accuracy, which is for example realized 
by the 10 12 of 3 He 2+ + e~ — >• T + + v e scattering events or determining the decay width of another nuclear species 
with 10 -6 accuracy, then, we can determine IV^I with 10 -6 accuracy. Alternatively, if we have the ability to detect 
the 0.02 eV m„, we can determine Ca with 10~ 6 accuracy. Then we can determine \V u d\ with 10~ 6 accuracy. 

According to Ref. 0, today's experimental value is \V u d\ — 0.97418 ± 0.00027. Hence, this method may have a 
great impact not only on the lepton sector but also on the quark sector. 



V. SUMMARY AND DISCUSSION 



We explained the new method to detect the neutrino mass. This method asks the precision measurement of mass 
defect and the released thermal energy. The results are written up in TABLE |TJ 

We also pointed out that this method is useful for \V u d\ precision measurement. This is because we can determine 
r precisely in this method. 

In TABLE HI the error of mx — "t He + — m e will reduce to 10 -6 in near future as explained in Appendix [Dj Also, 
the error of Mq can be already reduced to 10 -8 by the present technology (Appendix IB1 . 

In Eq. (|4]) and then also in Eqs. ([5]) and (fT9)) , we dealt with Mi — Mf. It is modified by considering the heat 
capacity of the sample. However, its effect can be controlled easily if the temperature of the sample is kept in low as 
explained in Appendix [E] 

Today, the experiments of elementary particle physics are categorized as the accelerator or non- accelerator ones. In 
both of them, the detected quantities are energy and momentum of the particle, which are the microscopic quantities. 
However, our method is an experiment which reveals the properties of elementary particles using the macroscopic 
quantities. 

The normalized standard deviation of (E v ) in one beta-decay event is o / (E v ) — 0.303 and for 1 mol tritium beta 
decays, it becomes 



(20) 
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where Na = 6.02 x 10 23 [mol x ] is the Avogadro constant. Therefore, we do not have to consider seriously the 
statistical error for 0.2 eV neutrino mass. Our work is unique in this respect. 
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Appendix A: Effective Neutrino Mass 



The decay width written in the mass eigenstates nik is approximated as 



3 3 3 

r = J2 \Uek\ 2 T(m k ) ~ J2 IM 2r (°) + H We. 



m 2 k d 2 T(m k ) 



ek] 2 dm 2 , |mfc=0 

k=l k=l k=l k 

rm\ , ml d 2 T(m v ) 

r (°) + — dm 2 U*=o - r(m„), 



(Al) 



where U e k are the leptonic mixing matrix elements, and 

3 3 

J2\U ek \ 2 = l, Y,\ U ^ 2m k^ m l- ( A2 ) 

k=l k=l 

As you can see from Eqs. ([6]) and ©, the linear term does not appear. 

Appendix B: Heat Value of the Source 

When a T decays into 3 He + + e~ + D e , the electron and 3 He + have about 6 keV momentum in total. This becomes 
the thermal energy in source. If we suppose that the pure tritium source has 3 g mass, it has about =6.02 x 10 23 
tritium atoms. According to the tritium half-life (12.33 years), the heat value in one second is 

N (l - e- r *)| 4=1[scc] x 6[keV] ~ 1[J]. (Bl) 
Ref. Q explains that the differential scanning calorimetry has the 25 [nW] resolution. This is 0(1O -8 ) of 1 [J/s]. 



Appendix C: Suppression of Ca Dependence in (E v ) 

Here, we show how the Ca dependence in (E v ) is suppressed. For simplicity, we represent the differential decay 
width as 

*L „ {aEl + b^) + C A (cEl + d-i), (CI) 
dhj v rriT Wt 
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where the coefficients a, b, c, and d are in the same order. Hence, the integral is easy to evaluate: 

dE ® ~ (^ Ma * 3 + \^—) + C A ( C -E^ 3 ' d E » ^ 



v 3 " 4 m T 3 " 4 m T 

Max 



5 hit 

Then, for E^f ax /mT <C 1, the mean neutrino energy becomes 



In tritium beta decay, the suppression factor E^ Lax /(20rrn) is O(10 6 ). 



(C2) 



JdEyE^ 3 b + dC A l 

^ - -jie^t 1 {1+ ^cu2o^ ] - (C3) 



Appendix D: mr — m Hc + — m e 

Ref. Q says that "In the history of mass spectrometry the precision of atomic mass determination has shown a 
constant improvement of about an order of magnitude every decade" . Also ref. fiol ] reported that the mass deference 
between tritium and helium-3 is 18590.1(1.7) eV in 1993. These facts suggest that we will be able to determine 
tot — m H_ c + ~ m e with 10 -6 error in 2013. Actually, the Max Planck Institute for Nuclear Physics aims for 10~ 7 
according to Ref. [TTj | . 



Appendix E: Heat Capacity 

If we consider the heat capacity of the sample, the equation Q is replaced by 

Ml-M' 

N= — (El) 

TOT - TORe 

where M[ and Mj: are the masses of initial and final state sample which are defined at zero temperature, respectively. 
These have the relations, Mi — M[ + C{T and Mf = Mj + CiT, where T is the sample temperature, and C and C/ 
are the sample heat capacity in the initial and final state, respectively. Then, the mean neutrino energy is given by 

lw\ ( \ft M a-TAC\ 

{E u ) = (m T - TO Ho ) 1 

(E2) 



AM - TAC J 

where AC = C\ — Cf. In the second line, we have used an approximation, TAC/ 'AM <C 1. 
Also, the equation (fT9l becomes 

= M[ - M< 2 ^ AM 2 f TAC, TAC,\ 

M[ - M' fl AMi ^ AA/i AM 2 J' { ' 

where the indices 1 and 2 represent the values at time t\ and t 2 , respectively. 

We can estimate the magnitude of these correction terms according to the Debye model. The Debye temperature 
Td of H2 is 105 K. Then, in low temperature (T <C Tb), the specific heat is expressed as 

C v 234 x N A k B ( ^ 234 x (6.02 x 10 23 [mob 1 ]) x (1.38 x 10~ 23 [J/K]) ( (7^-7) 



— J ~ 234 X X 11T" [mol -jj X ^1.38 X 1U [J/Kjj ^ — J y — J 



- 0.168 x ^—j [J/(K-mol)], 

where k B — 1-38 x 10~ 23 [J/K] is the Boltzmann constant and N A — 6.02 x 10 23 [mol -1 ] is the Avogadro constant. 
5.5 % of tritium decay in one year since the half-life of tritium is about 12.32 year. Assuming that tritium and 
helium-3 have the same order of specific heat, we give AC ~ 0.168 [J/(K • mol)] x 1 [mol] x 0.055 ~ 0.01 [J/K]. Also, 
AM ~ 18keV x N A x 0.055 x (1.6 x 10~ 19 [J/eV]) ~ 10 8 J. Then TAC/ AM ~ 10[K] x 0.01[J/K]/10 8 [J] ~ 10~ 9 . This 
means that we can easily estimate the heat capacity effect enough accurately in the low temperature system. 



